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An operator-theoretic approach to differential positivity 

A. Mauroy, F. Forni, and R. Sepulchre 


Abstract —Differentially positive systems are systems whose 
linearization along trajectories is positive. Under mild as- 
snmptions, their solutions asymptotically converge to a one¬ 
dimensional attractor, which must be a limit cycle in the absence 
of fixed points in the limit set. In this paper, we investigate 
the general connections between the (geometric) properties 
of differentially positive systems and the (spectral) properties 
of the Koopman operator. In particular, we obtain converse 
results for differential positivity, showing for instance that any 
hyperbolic limit cycle is differentially positive in its basin of 
attraction. We also provide the construction of a contracting 
cone field. 

I. Introduction 

A linear system is positive if its trajectories leave some 
conic subset JC of the state space invariant [3]. Positivity is 
at the core of a number of applications because it strongly 
restricts the linear behavior [5], [16], [17]. Under mild condi¬ 
tions, Perron-Frobenius theory guarantees that every bounded 
trajectory converges asymptotically to a one-dimensional 
attractor given by the ray Awe/C, where AgK and w is 
the dominant eigenvector of the system state matrix [2], [3]. 

Differential positivity brings linear positivity to the nonlin¬ 
ear setting. A nonlinear system is differentially positive if its 
linearization along any trajectory leaves some cone (held) 
invariant [8]. Differentially positive systems are a sizable 
class of systems, which includes monotone systems [1], 
[9], [19]. Differential positivity also restricts the asymptotic 
behavior of a nonlinear system. Under mild conditions, a 
suitable differential formulation of Perron-Frobenius theory 
guarantees that the trajectories of the nonlinear systems 
converge asymptotically to a one-dimensional attractor. In 
contrast to linear positivity, this attractor is not necessarily a 
ray but a curve, possibly given by a collection of hxed points 
and connecting arcs or by a limit cycle [8], [7]. 

Differential positivity is a promising tool for the study of 
bistable and periodic behaviors, since it reduces the analysis 
of those behaviors to the characterization of a suitable cone 
held on the system state manifold. However, besides spe- 
cihc families of monotone systems, there is no constructive 
methodology to hnd such cone helds. At more fundamental 
level, it is not even clear how demanding is to use differential 
positivity for capturing bistable and periodic behaviors. This 
paper provides a hrst answer to both these questions. 

A. Mauroy is with the Department of Electrical Engineering 
and Computer Science, University of Liege, 4000 Liege, Belgium, 
a.mauroy@ulg.ac,be 

F. Forni and R. Sepulchre are with the University of Cambridge, 
Department of Engineering, Trumpington Street, Cambridge CB2 IPZ, 
and with the Department of Electrical Engineering and Computer Sci¬ 
ence, University of Liege, 4000 Liege, Belgium, f f 2 8 6 @ cam. ac . uk, 
r.sepulchre@eng,cam.ac.uk 


Bridging the (geometric) properties of differentially pos¬ 
itive systems and the (spectral) properties of the Koopman 
operator [14], the paper illustrates the tight connection be¬ 
tween the existence of a suitable collection of Koopman 
eigenfunctions for the system and the construction of a 
cone held. This approach leads directly to a numerical tool 
for constructing cone helds. At more fundamental level, 
Koopman theory provides a way to derive converse results 
for differential positivity. The striking outcome is that any 
system with a hyperbolic limit cycle is differentially positive 
in the basin of attraction of the limit cycle. 

A short introduction to both differential positivity and 
Koopman operator theory is provided in the next section, 
which follows a brief discussion on the basic geometric 
tools used in the paper. The connection between Koopman 
operator theory and differential positivity is developed in 
Section We provide an explicit construction of the cone 
held based on a suitable set of Koopman eigenfunctions, 
and we show the precise relation between the so-called 
Perron-Frobenius vector held and the dominant Koopman 
eigenfunction. Section |iy] is dedicated to converse results 
for hyperbolic hxed points and hyperbolic limit cycles. 
Section [V] provides a cone held for a system with a stable 
equilibrium and for the Van der Pol oscillator, by exploiting 
numerical methods for computing Koopman eigenfunctions 
based on Laplace averages. Conclusions follow. Proofs are 
in appendix. 

II. A GLIMPSE INTO DIFFERENTIAL POSITIVITY AND 

Koopman operator 
A. Manifolds and prolonged dynamics 

The exposition of the paper takes advantage of a few basic 
geometric notions on Riemannian manifolds. Let A be a 
smooth n-dimensional manifold endowed with a Riemannian 
metric (•,-)x : T^X x T^X M where T^X denotes 
the tangent space at x G X. We will use |(5a;| to denote 
y {Sx, 5x)x for all 5x G T^X. Given two smooth manifolds 
Xi and X 2 , and a differentiable function (or observable) g : 
Xi -G X2, let dg{x) : T^Xi -G Tg^^^X2, 6x i-A dg{x)[Sx], 
be the differential of g at x. When clear from the context, 
we will simply write dg{x)6x = dg{x)[6x], and we will use 
dg(x)A := {dg{x)Sx\Sx G A} for all A C T^X. 

The paper focuses on continuous-time dynamical systems 
S on A represented by i = f{x), where x G X and f{x) G 
TxX. We assume that / G C'^(A) (twice differentiable) and 
that the system is forward and backward complete, that is, the 
How '0:KxA—J-AofE satishes = f{ip(t,x)) 

for alH € K and a; S A. In what follows, for simplicity, we 
will also use the mapping '0‘(-) = '0(f, •) : A —> A and we 



will sometimes refer to the trajectory x{-) = '(/>(•, a:o) of S 
from the initial condition x{to) = tjj{tQ,xo). 

To characterize the property of differential positivity we 
will make use of the notion of prolonged dynamics of 
S, represented by 

Ay • / ^ ^ 

' \ Ax = df{x)Sx 

where (x, Ax) belongs to the tangent bundle TX = 
[JxGxi^} ^ T'xX, [4]. The flow of AS is a mapping in 
K X TX TX specified by (x,Ax) !->■ {ip^{x),d'4>^{x)5x). 


B. Differential positivity 

A linear system on M" is positive if there exists a cone 
JC C M" which is forward invariant for the system dynamics. 
Indeed, given x = Ax, positivity reads C /C for all t > 
0. Differential positivity is a way to extend linear positivity 
to nonlinear dynamics, by requiring that a given cone (field) 
is forward invariant for the prolonged dynamics AS, [8]. 

We endow the manifold X with a cone field 

/C(x) C TxX Vx G A . 

Each cone K.{x) is closed and solid, and satisfies the follow¬ 
ing properties: for all x € X, (i) /C(x) + /C(x) C /C(x), (ii) 
alC{x) C /C(x) for all a G K'*", (iii) /C(x) D —/C(x) = {0} 
(i.e. convex and pointed). The forward invariance of the cone 
field along the prolonged dynamics AS reads as follows [8]. 

Definition 1 (Differential positivity): The system S is dif¬ 
ferentially positive (with respect to the cone field /C) if the 
flow of the prolonged system AS leaves the cone invariant 

dip*(yx)JC{x) C /C(^*(x)) Vx G <T , Vf > 0. 

In addition, S is (uniformly) strictly differentially positive if 
it is differentially positive and if there exist a constant T > 0 
and a cone field TZ{x) C int/C(x) U {0} such that 

d'ip^(x)/C(x) C 7Z(tp^(x)) Vx e X ,Vt > T. ^ 

An illustration of the strict differential positivity property is 
provided in Figure [T] 
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Fig. I. Differential positivity requires the forward invariance of the cone 
field IC{x) along trajectories. Strict differential positivity requires that along 
trajectories the rays of the initial cone converge towards each other. 


To avoid pathological cases, we assume that for every pair 
of points xi, X 2 G X, there exists a linear invertible mapping 
r(xi,X 2 ) : Tx^X —)• Tx^X such that r(xi,X 2 )/C(xi) = 


JC{x 2 ) and r(xi,X2)7^(xi) = 72.(x2). Furthermore, we 
consider the representation 


lC{x) = {Sx G TxX \ ki{x,Sx) > 0 Vi = l,...,m} 

TZix) = {Ax G I fci(x, ilfy) > e Vi=l,...,TO} 

for some m G N and e > 0, where : FA —>■ K are smooth 
functions (the reader is referred to [8], [7] for details). 

It follows from Definition[T]that E is differentially positive 
if from any initial condition ki{x, Ax) > 0, z G (1, ..., m}, 
the prolonged system satisfies ki{4’^{x),dff{x)5x) > 0 for 
all i G to} and all t > 0. In addition, strictly 

differential positivity requires [^\x), > e, for 

all f > r. In compact sets, both properties can be checked 
by simple geometric conditions [7]. 

Strictly differentially positive systems enjoy a projective 
contraction property [8], [3], which leads to the existence of 
the so-called Perron-Frobenius vector field w(x) G int/C(x), 
the differential equivalent of the Perron-Frobenius eigenvec¬ 
tor of linear positive mappings. The ray {Aw(x) | A > 0} C 
int/C(x) U {0} is an attractor for the prolonged dynamics, in 
the precise sense that 


dip\ip *(x))Ax ^ 

i-voo 


( 1 ) 


for all X G A and Ax G lC(ff~*(^x)) \ {0}. It follows that, 
w(z/^*(x)) = d'ip^{x)'w{x)/\d'ip*{x)'w{x)\. 

In the next sections we will use the notion of Perron- 
Frobenius curve : / C K —A, which is an integral 
curve of the Perron Frobenius vector field, i.e. = 

w( 7 '^(s)) for all s G /. 


C. Koopman operator 

The so-called Koopman operator describes the evolution 
of observables g ■. X ^ along the trajectories of E. 

Definition 2 (Koopman operator): For a given functional 
space Q, the (semi-)group of Koopman operators U* : Q ^ Q 
associated with a system E is defined by 

U*g = goijjff g gG ,tGR 

where V’*(') is the flow map of E. o 

Even when E is a nonlinear system, the Koopman operator 
is linear, so that it can be studied through its spectral 
properties. 

Definition 3 (Koopman eigenfunction and eigenvalue): 
The observable ()>a G is an eigenfunction of the Koopman 
operator (called Koopman eigenfunction hereafter) if there 
exists a value A G C such that 


{(pxofi^fix) = U*-4)\(x) = e^^(l)x{x) 'ixGX. (2) 


The value A is the associated eigenvalue. o 

Remark 1: When G P= C'^(A’), the semi-group of Koop¬ 
man operators admits the infinitesimal generator L : G G 
defined by Lg{x) = dg{x)[f{x)] and we have U* = In 
this case, a Koopman eigenfunction fix satisfies 

Lfixix) = d(j)x{x)[f{x)\ = Mxix) ■ (3) 

o 





The Koopman eigenfunctions capture important geometric 
properties of the dynamics (see e.g. [12], [13], [14]). In 
Section III they will be related to the differential positivity 
properties of the system. As a preliminary, we consider the 
semi-group of Koopman operators C/* associated with the 
prolonged system (5S, i.e. 

U*g{x,Sx) = g{'ip*{x),d'ip*{x)6x) 


with the observables g : TX —> C. We have the following 
result. 

Lemma 1: Suppose that 4)\ G C^(X) is an eigenfunction 
of {7* associated with the system E. Then the Koopman 
operator U* associated with the prolonged system 5E admits 
the eigenfunctions (in the appropriate functional space) 

^^^\x,5x) = (j}x{x) 

^^^\x,5x) = d(j)\(x)5x 

for all {x, 6x) gTX. o 

Proof: It is clear that In addition, 

U*f^^\x,5x) = {x))[d'il)*' {x)5x\ 

= d{(j)xO'f*){x)5x 
= e^^d(j)\{x)6x 
= e^*^^x\x,6x) , 

where we used 0 . ■ 

If E admits a fixed point, the prolonged system i5E is charac¬ 
terized by a star node (with eigenvalues of multiplicity two). 
In this case, it is known that the corresponding Koopman 
eigenvalues are also of multiplicity two (see Remark 2 in 
[13]). 


III. From Koopman eigenfunctions to cone eields 

In this section, we present general results that connect dif¬ 
ferential positivity to the spectral properties of the Koopman 
operator. We show that a system is differentially positive if 
there exist specific independent Koopman eigenfunctions. 

Proposition 1: Suppose that the n-dimensional system E 
admits a set of Koopman eigenfunctions fx^ G C^{X), j = 
1 ,..., n, [ft{Ai} > 3fi{Aj}, such that the linear map : 

T^X -)■ C", 

i9<I>(a;)^a; = {dfxi ix)Sx ,..., ix)dx) 


Every cone field JC is the local representation of a (global) 
conal order ^ CXxX. The order ^ is derived from JC by 
integration: xi -< X 2 if and only if there exists a curve 7 : 
[si, S 2 ] C M —)• A" with 7(si) = xi and 7(52) = X 2 such that 
2^7(3) G /C(j(s)) for all s G [si, S2]. It is noticeable that the 
conal order given by the cone field of Proposition [1] has the 
following equivalent characterization: xi^X 2 if and only if 


^Ai(a: 2 ) - (/>\i(xi) + \(t>\j{x 2 ) - (t)Xj{xi)\>0, for all j>2 
(replace fxi by in case 2)). Indeed, the contracting 

cone field defined locally by the Koopman eigenfunctions 
of the prolonged system (5E induces a conal order which is 
captured directly by the Koopman eigenfunctions of E. 

Remaric 2: When 3?{Aj} < 0 for all j in Proposition 
the Koopman eigenfunctions of 5E can be used to construct 
a differential Finsler-Lyapunov function which decays along 
the trajectories of the system. By integration, a differential 
Finsler-Lyapunov function induces a distance on the system 
state manifold, which also decays along any pair of tra¬ 
jectories of E, establishing contraction [ 6 ]. As above, it is 
noticeable that such distance has an equivalent characteri¬ 
zation based on the Koopman eigenfunctions of E. Indeed, 
the Koopman eigenfunctions of E capture the contractive 
behavior of the system [11], [13] 

The next proposition illustrates the relationship between 
Koopman eigenfunctions and the Perron-Frobenius vector 
field. 

Proposition 2: Suppose that there exists a set of eigen¬ 
functions (pXj that satisfies the conditions of Proposition ^ 
with 5ft{Ai} > 5ft{Aj} for all j (i.e. E is strictly differentially 
positive). Then the Perron-Frobenius vector field is the 
unique vector field w : X ^ T^X, |w(-)| = 1, that satisfies 


d(j)Xj{x)[w{x)] = 0 j >2. (7) 

o 

Proposition implies that the Perron-Frobenius vector 
field is related to zero level sets of Koopman eigenfunctions 
associated with the prolonged system i5E. At global level, 
the integral curves of the Perron-Frobenius vector field 
correspond to the intersection of the level sets of Koopman 
eigenfunctions associated with E. Precisely, any set 

m 

{^{xGX\fxM = C,} , (C'2,...,C™)gC"-i (8) 

i=2 

is the image of some Perron-Frobenius curve. 


is injective for all x G X. Then E is differentially positive 
if one of the following is satisfied: 

1) Ai G K. The cone field reads JC{x) 

{SxGTxX\ dfxAx)^x - \d(t)Xj{x)5x\ > 0, Vj > 2} ; 

(5) 

2) Ai G with Zfxi G C^{X) and I^aJ is constant on 
X. The cone field reads IC{x) 

{Sx G TxX I dZfxi {x)Sx - \d(j)Xj (a::)(5a:| > 0, Vj > 2} . 

(6) 

The system is strictly differentially positive if 5ft{Ai}> 3fi{A2}. 


IV. Converse results eor dieferential positivity 

Restricting the analysis to systems E on vector spaces 
X := R”, we show that the presence of a stable hyperbolic 
fixed point or of a stable hyperbolic limit cycle is a sufficient 
condition for E to be strictly differentially positive in their 
basin of attraction. 

Proposition 3: [Hyperbolic stable fixed point] Consider a 
system x = f{x), with x G R" and / G C^, which admits a 
fixed point x* with a basin of attraction B(x*) C R". For 
j = 1,... ,n, assume that the eigenvalues Xj of the Jacobian 
matrix J = § 7 ( 2 :*) satisfy 0 > 5ft{Aj} > JilAj+i} and that 
the eigenvectors are independent. The system is differentially 




positive in B{x*) if and only if Ai S M. Moreover, it is 
strictly differentially positive if Ai > 5ft{A2}. o 

The result of Proposition also holds if the fixed point 
is unstable with 3?{Aj} > 0 for all j but it does not 
hold with a saddle node. This is not surprising since [8, 
Corollary 3] shows examples of hyperbolic saddle nodes that 
are incompatible with differential positivity. 

Proposition 4: [Hyperbolic stable limit cycle] If a system 
X = f{x), with a; S K" and / S C^, admits a stable 
hyperbolic limit cycle T (with independent eigenvectors 
of the monodromy matrix), then it is strictly differentially 
positive in the basin of attraction yB(r) C M" of T. o 

Remark 3: Proposition and ([^ show that the Perron- 
Frobenius curves are the intersections of the level sets of n— 1 
Koopman eigenfunctions j = 2,... ,n. Therefore, the 
Perron-Frobenius curves can be interpreted as dual quantities 
to the (n — 1)-dimensional level sets of the Koopman 
eigenfunction the so-called isostables for fixed points 
[13] and the isochrons for limit cycles [12]. o 

V. Numerical computation 


B. Fixed points 

When the system admits a stable hyperbolic fixed point 
X*, a contracting cone field is given by <0, where the 
Koopman eigenfunctions are given by the Laplace averages 
Figure shows the contracting cone field for the 
dynamics 

xi = -sin(xi)-I-cos(x 2 ) - 1 
±2 = — cos(xi) — 1.5sin(x2) + 1 

which has a stable fixed point at the origin. 

The Laplace averages and (where Ai = — 1 and 
A 2 = —1.5 are the eigenvalues of the Jacobian matrix J at 
the origin) are computed with the observables g{x,6x) = 
vf6x and g{x,6x) = v^Sx, respectively, where vx and V 2 
are the left eigenvectors of J. This choice ensures that each 
average is finite (provided also that 2Ai < A 2 < Ai, 
which is the case here) and nonzero. Note also that we have 
g = dg with g = vj x, as required by Lemma 0 The Perron- 
Frobenius vector field is given by 0. An illustration is in 
Figure 


A. Preliminaries on Laplace averages 

In this section we exploit the theoretical results of the 
paper to derive contracting cone fields for fixed points and 
limit cycles based on the Koopman eigenfunctions of the 
prolonged system. 

When the trajectories of the system are available, an 
efficient method for computing Koopman eigenfunctions is 
based on Laplace averages (see e.g. [15]) 

gT{x)=lua I- f [go'iP%x)e~^^ dt ( 9 ) 

T-i-oo T Jq 

for an observable g : X ^ C. The average 0 is well- 
defined (i.e. finite) when A is a Koopman eigenvalue and g 
is a well-chosen observable (for instance, g must be zero on 
the attractor when 3?{A} < 0). From 0 , it is easy to see 
that U*g'^ = Provided that p™ 7 ^ 0, we can define 

= 9T- 

We can also obtain the Koopman eigenfunction 
^^^\x,Sx) = d4>\{x)Sx associated with the prolonged 
system ST, by computing the Laplace averages along the 
trajectories of i5S. This is summarized in the following 
lemma. 

Lemma 2: Suppose that the Laplace average is 

finite and nonzero, so that (j>\ = g^^. Then, we have 

^x\x,Sx) = gl'’{x,5x) , 

with 

1 

57 {x,fe)= lim - / g{ip\x),d'ip\x)5x)e dt 
T->00 T Jq 

( 10 ) 

and with the observable g{x,Sx) = dg{x)5x. o 

Proof: The result is obtained by differentiating 0. ■ 
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Fig. 2. Cone field (red) and Pen'on Frobenius vector field (blue) for GD- 


C. Limit cycles 

When the system admits a stable hyperbolic limit cycle, a 
contracting cone field is defined by 0 and can be expressed 
in terms of Koopman eigenfunctions associated with the sys¬ 
tem T and the prolonged system ST. The Laplace averages 
0 and ( [Tol l can be used to compute the cone field. 

Figure [^ shows the contracting cone field for the Van der 
Pol dynamics 

^ 2 ^ (12) 

X2 = {l-xf)X2-Xi 

which has a stable limit cycle F (of period T > 0). 

The Laplace averages p™ and g'Jf (with Ai = VI-k jT') 
are computed with the observables g{x) — [1 0]^x and 
g{x,Sx) — dg{x)Sx = [1 Oj^fe, respectively. According 
to ([T^, we have dTfxi = 'Tl'® average g'Jf 

(where A 2 is the nonzero Floquet exponent of the limit cycle) 
is computed with the observable g{x,Sx) = ^{p{x))'^Sx, 
where ^ ; F —> TA” is a unit vector field such that ^{x) is 
perpendicular to the tangent direction to the limit cycle at x, 
and where p : A” —F is a radial projection on F (i.e. p{x) 









is the intersection between T and the line passing through 
X and the origin). Note that, as required by Lemma HI g is 
the differential of an observable measuring a distance to the 
limit cycle. The Perron-Frobenius vector field is given by 
(0. An illustration is in Figure 
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Fig. 3. Cone field (red) and Perron Frobenius vector field (blue) for (ID- 

Remark 4: Besides the Laplace averages, Koopman eigen¬ 
functions can be computed by other methods, which also 
provide novel ways to derive cone fields. For example, the 
Koopman operator can be expanded on a (finite) polynomial 
basis [11], yielding a polynomial approximation for eigen¬ 
functions and cone fields. In a similar way, the extended 
dynamic mode decomposition [18] could also be employed. 

VI. Conclusion 


\d(j)\j{x)5x\ — \d(j)\.{x)5x’\ > 0 if Sx,Sx' S JC{x). 
(ii) alC{x) C lC{x), a > 0 since d(j)\^{x)[a5x\ — 
\d(t>\j{x)[a5x\\ > 0 if (5a; € /C(x). (iii) JC{x)f^—JC{x) = {0} 
since d(j)\^{x)5x — \d(j)\.{x)5x\ > 0 and —d(j)\i{x)Sx — 
\d(j)\.{x)6x\ > 0 imply d(j)\.{x)dx = 0 Vj, so that 6x = 0 
since the linear map i9<I>(a;) is injective. 

Using the Koopman eigenfunctions of the prolonged sys¬ 
tem (5S and Lemma [T] we have 

d(j)Xi{ip\x),dip\x)5x) - \d(j)x.{il)\x),di)\x)5x)\ 

= U^^^^^ix,Sx) - \ij^^^^^{x,6x)\ 

= Sx) — i5x)|^ (13) 

> (d(/>Xi(x)Sx - ldc/^Xj(x)Sxl) 

> 0 

for all a; € A” and Sx G K.(x) and for all t > 0. It follows 
that d'ilj\x)JC{x) C /C(V'‘(a:)). If 3fJ{Ai} > 5R{A2}, ([B]) 
is a strict inequality and E is strictly differentially positive 
(uniformly with an arbitrary T > 0). 

Case Ai G iM. We consider the cone field (|^. It follows on 
similar lines that (i) IC{x) +JC{x) C ]C{x) and (ii) alC{x) C 
IC{x), a > 0. In addition, dZ(j)Xi{x)6x — \d<j)Xj{x)6x\ > 0 
and —dZ(j)Xi {x)6x—\d(j)Xj (a;)(5a;| > 0 imply dZ(j)Xi {x)6x = 
d(j)Xj{x)5x = 0 Vj > 2. Since is constant by 

assumption, we have 

d(t>Xi = 9(|(;!)Ai|e*^'^^i) = (14) 


Reducing the analysis of a nonlinear finite-dimensional 
system to the analysis of an infinite-dimensional linear sys¬ 
tem, Koopman operator theory is a powerful tool of nonlinear 
control. In this paper we bridged Koopman operator theory 
and differential positivity, opening the way to the use of 
spectral methods for differential positivity. We illustrated 
a tight relation between Koopman eigenfunctions and cone 
fields, leading to converse results for differential positivity 
and to numerical tools for the construction of cone fields. 

The bridge between Koopman theory and differential 
positivity opens new interesting directions of research. For 
example, we believe that the approach pursued in the paper 
will lead to converse results for the larger class of normally 
hyperbolic one-dimensional attractors. This will be the object 
of future research. 
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Appendix 

Proof of Proposition 

Case Ai G K. We first show that 0 is a well-defined 
cone field, (i) /C(x) + IC{x) C JC{x) since d(pxiix)[dx + 
5x'] — \d(j)Xj{x)[5x + (5x']| > d(j)Xiix)Sx + dfxi^Sx' — 


and it follows that d(j)Xi{x)Sx = 0. Then the injectivity of 
i9<l>(x) implies (5x = 0 , so that IC{x) n —/C(x) = {0}. 

Since \(j)Xi \ must be nonzero, it follows from ([T^ that 


d(l)Xrix)Sx _ fx^ix,6x) 


dZfxi (x)Sx = 


Using Lemma and 3?{Aj} < [ft{Ai} = 0, we have 

5Z(/)Ai(i/’‘(a;))[5i/)*(x)(5x] - \d(l>Xj{'f\x))[df*{x)6x]\ 
_ (x, (5x) 

iWfxiix) 

_ _ SR{A,}i| T(2)/ ^ ^1 

> dZfxi {x)6x - \d(j)Xj (x)fe| 

> 0 


(15) 


- \U^f^^HxAx)\ 


for all X G Y and 5x G /C(x) and for all t > 0. It follows that 
dtp*{x)IC{x) C K.{tjf{x)). The proof for strictly differential 
positivity follows on similar lines. ■ 

Proof of Proposition 

Consider o; G K such that jO^Ai (x)[aw(x)]| = 1. If Ai G 
K, this implies that c)(^Ai(x)[q;w(x)] has a unique value up 
to a sign. If Ai G iR, |()>AiI is constant and ( [T^ implies that 
i9(^Ai(x)[aw(x)] has a unique value up to a sign. Since ST* 
is injective, it follows that w(x) is unique (up to a sign). 







Next, we show that the Perron-Frobenius vector field 0 
satisfies 0. We have 

\d<i)x^ {x)[dil}\'tp~\x))8x\\ 

= \U^^x^,i'>l’~\x),5x)\ ( 16 ) 

where we used x = ip*i^)) and Lemma If Ai € K 
and Sx G IC{ip~*{x)), it follows from ( [T6] l that 

|i9<^A,(a;)[i9V'‘(V’"‘(a;))i5a;]| 

< e^^^^^*d(Px^i'tp~\x))Sx 

= i^-\x),Sx) 

= e^^^^^^-^^^*^(|)x^ix)[^'lP\PJ-\x))6x] 

< (^)ll \dpj\ip-\x))Sx\ 

with j > 2 and with ||9(/);^j(a;)|| = max| 53 .|^i \d(j)Xi{x)Sx\. 
In the case Ai G jM, it follows on similar lines that 

\d(px, {x)[d'4!\tp~\x))Sx]\ 


Assuming without loss of generality that {(pxiix)] = 1 when 
Ai G iM, ( pAl i implies that ||9Z());^J| = ||9());vJ|. Then, in 
both cases Ai G K and Ai G iK, we have 


\d<px,{x)[d'ip*{pj *(a^))fa]| 
\dtp*{'ip~*{x))6x\ 




Finally, taking the limit t ^ oo and using the definition 0, 
we obtain (a;)[w(x)] | = 0 since 3?{Aj} < 3?{Ai}. This 
concludes the proof. ■ 


Proof of Proposition 

Sufficiency. It is known that the system admits n indepen¬ 
dent Koopman eigenfunctions fxj € C^{B{x*y) associated 
with the eigenvalues Aj of J. This result follows from 
Theorem 2.3 in [10], which shows the existence of a 
diffeomorphism h : B{x*) —>■ K” such that y = h{x) and 
y = Jy. Letting fxjix) = vjh{x) where vj is the left 
eigenvector of J associated with Xj, we verify that 

d(px,{x)[f{x)] = vjdh{x)[f{x)] = vjjy = X^vjy = Xjfx, 

and 0 implies that fxj is a Koopman eigenfunction. In 
addition, since is a diffeomorphism, dh is injective and the 
linear map d^{x) = {vfdh{x),... ,v^dh{x)) is injective 
since the eigenvectors Vj are independent. Then the result 
follows from Propositionwith X = B{x*) and Ai G K. 

Necessity. If the system is (strictly) differentially positive, 
the differential dynamics 5x = J 5x at the fixed point 
is (strictly) differentially positive. Then the invariant cone 
IC{x*) must contain the dominant direction, which corre¬ 
sponds to the right eigenvector v of J associated with Ai. 
If Ai f. K, it is clear that the invariance of the cone 


field implies that JC{x*) contains the entire two-dimensional 
plane spanned by [ftlf;} and S5{u}. Thus we have IC{x*) n 
—K.{x*) f {0} and /C(a;*) is not a pointed cone. This is a 
contradiction, so that Ai G M. ■ 


Proof of Proposition 

It is known that the system admits n independent Koop¬ 
man eigenfunctions fxj G C^{B{x*)). This result follows 
from Theorem 2.6 in [10] (and from Floquet theory), which 
shows the existence of a diffeomorphism h = {hy, he) : 
B{T) -G X such that {y,0) = {hy{x),he{x)) 

and y = By, 9 = (jj, where uj = 2 t:/T (with T the 
period of the limit cycle). It is clear that fx^ix) = 
is a Koopman eigenfunction associated with the eigenvalue 
Ai = iuj. In addition, the matrix B has n — 1 eigenvalues 
Xj, j = 2,...,n (i.e. the nonzero Floquet exponents of 
the limit cycle) associated with the left eigenvectors vj. 
Since the limit cycle is hyperbolic and stable, we have 
3fJ{Aj} < 0 for all j > 2. As in the proof of Proposition 
^ it follows that fxjix) = vjhy(x) are Koopman eigen¬ 
functions associated with the eigenvalues Xj. In addition, 
since is a diffeomorphism, dh is injective and the linear 
map d^{x) = {ie^^'^''^^dhe{x),V 2 dhy{x),... ,v^dhy(x)) 
is injective since the eigenvectors vj are independent. Finally, 
Z(/)Ai = he G and I^aJ = 1- Then the result follows 
from Propositionwith X = B{x*) and Ai G jM. ■ 
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